#pragma rtGlobals=1		// Use modern global access method.

macro PLhosei()
	PLhosei1()
end
function PLhosei1()

	string xdata
	prompt xdata,"Choose X(nm) data.",popup,wavelist("*",";","")
	string ydata
	prompt ydata,"Choose Y data.", popup, wavelist("*",";","")
	doprompt "execute",xdata,ydata
	
	if(V_Flag!=1)
		if(strlen(ydata)==0 || strlen(xdata)==0)
			doalert 0,"it failed !"
		else
			PLchange(xdata,ydata)
		endif
	endif
end

function PLchange(xdata,ydata)
	
	string xdata
	string ydata
	variable i=0
	wave xnm=$xdata
	wave yvalue=$ydata	
	
	duplicate xnm, xev, deltaEV, deltaNM
	duplicate yvalue, yvalueEV
	
	//change nm to ev for X data
	do
		xev[i]=1240/xnm[i]
		i=i+1
	while(i<numpnts(xev))
	i=0
	//make delta X value
	do 
		if (i==numpnts(xev)-1)
			deltaEV[i]=abs(xev[i]-xev[i-1])
			deltaNM[i]=abs(xnm[i]-xnm[i-1])
		else
			deltaEV[i]=abs(xev[i+1]-xev[i])
			deltaNM[i]=abs(xnm[i+1]-xnm[i])
		endif
		i=i+1
	while(i<numpnts(xev))
	i=0
	//Change yvalue(nm) to yvalue(ev)
	do
		yvalueEV[i]=(yvalue[i]*deltaNM[i])/deltaEV[i]
		yvalueEV[i]=yvalueEV[i]/1000000
		i=i+1
	while(i<numpnts(yvalue))
	
	duplicate/o xev,  $xdata + "_ev"
	duplicate/o yvalueEV, $ydata + "_EV"
	display $ydata+"_EV" vs $xdata+"_ev"
	killwaves/z deltaEV, deltaNM, xev, yvalueEV
		
End

//from here, chokusenfitfunction
macro CHOKUSENdeHIKU()
	chokusenhiku1()
end

function chokusenhiku1()

	variable i=0, s=0, t=0
	
	if(strlen(csrwave(a)) == 0 || strlen(csrwave(b))==0)
		doalert 0,"threre is no cursor."
	else
		duplicate/o csrwaveref(a), tempw
		duplicate/o csrwaveref(a), tempw2
		s=(vcsr(a)-vcsr(b))/(hcsr(a)-hcsr(b))
		t=(hcsr(a)*vcsr(b)-hcsr(b)*vcsr(a))/(hcsr(a)-hcsr(b))
		do
			tempw[i]=s*csrxwaveref(a)[i]+t
			tempw2[i]=csrwaveref(a)[i]-tempw[i]
			i=i+1
		while(i<numpnts(csrwaveref(a)))
		
		duplicate/o tempw2, $csrwave(a)+"_2"
		appendtograph $csrwave(a)+"_2" vs csrxwaveref(a)
		killwaves/z tempw2, tempw
	endif
			
end

//from here, threre are integral macros
macro SEKIBUNsuru()
	sekibun1()
end

function sekibun1()
	
	variable i=0, s=0, t=0, menseki=0, j=0
	
	if(strlen(csrwave(a)) == 0 || strlen(csrwave(b))==0)
		doalert 0,"threre is no cursor."
	else 
		duplicate/o csrwaveref(a), tempw
		duplicate/o csrwaveref(a), tempw2
		s=(vcsr(a)-vcsr(b))/(hcsr(a)-hcsr(b))
		t=(hcsr(a)*vcsr(b)-hcsr(b)*vcsr(a))/(hcsr(a)-hcsr(b))
		do
			tempw[i]=s*csrxwaveref(a)[i]+t
			tempw2[i]=csrwaveref(a)[i]-tempw[i]
			i=i+1
		while(i<numpnts(csrwaveref(a)))
		
		if(pcsr(a) >= pcsr(b))
			j=pcsr(a)
			i=pcsr(b)
		else
			j=pcsr(b)
			i=pcsr(a)
		endif
		
		do
			menseki=menseki+csrwaveref(a)[i] - tempw[i]
			i=i+1
		while(i<j+1)
		
		print "+++++++++++++++++++++++++++"
		print  "sekibun value = " + num2str(menseki)
		print "+++++++++++++++++++++++++++"
		killwaves/z tempw2, tempw
		
	endif
	
end